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Abstract. We report the determination from high-resolution spectra of the atmospheric parameters and abundances of 13 
chemical species (among which lithium) in 8 single-lined active binaries. These data are combined with our previous results 
for 6 other RS CVn systems to examine a possible relationship between the photospheric abundance patterns and the stellar 
activity level. The stars analyzed are generally found to exhibit peculiar abundance ratios compared to inactive, galactic disk 
stars of similar metallicities. We argue that this behaviour is unlikely an artefact of errors in the determination of the atmospheric 
parameters or non-standard mixing processes along the red giant branch, but diagnoses instead the combined action of various 
physical processes related to activity. The most promising candidates are cool spot groups covering a very substantial fraction 
of the stellar photosphere or NLTE effects arising from nonthermal excitation. However, we cannot exclude the possibility 
that more general shortcomings in our understanding of K-type stars (e.g. inadequacies in the atmospheric models) also play a 
significant role. Lastly, we call attention to the unreliability of the (V - R) and (V - I) colour indices as temperature indicators 
in chromospherically active stars. 
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C3 1 ■ Introduction tent can also be used to set stringent constraints on the evolu- 

^ tionary status, but accurate effective temperatures are needed. 

£ ; The class of RS CVn binaries is loosely defined as being com- In yirtue of ^ growing body of evidence for distinct abun . 

^ . posed of systems made up of two late-type, chromospherically dance ^^ of SQme elements fa ^ and thick disk stars (e g 

in active stars, at least one of which having already evolved off the Bensbyi Fekzing & Lundstr5m 2003), elemental abundances 

.*? , main sequence (Hall 1976). Because the chromospheric activ- can alsQ complement the kinema tical data to help establishing 

ity is rejuvenated by the spinning up of the stars due to the whether these objects bdong tQ young Qr oM digk staf popuk _ 

transfer of orbital angular momentum to the rotational spin, t j nns 

they can sustain high activity levels at advanced ages com- „, , . _. , • . ., . • . , • • ., ■ 

3 fe 3 & We have started a project that aims at deriving the pho- 

pared to otherwise similar, single stars. A better understanding . , • . juj e iiui 

F ' & & tospheric parameters and abundances of several key chemical 

of this phenomenon, and in particular of the temporal evolu- ... , . . ,. . _. . „ TC 

F ' f t- species in a large sample of active binaries. First results for IS 

tion of the stellar activity level (e.g. the coronal X-ray emis- w and vg51 Cen haye been presented fa Katz et al (2003i 

sion), requires a good handle on their evolutionary status (e.g. hereafter paper r) Four additional RS CVn binaries have been 

Barrado et al. 1994). Precise estimates of the atmospheric pa- analyzed ^ ^ second paper of ^ geries (Mord ^ ^ 2QQ ^ 

rameters and abundances are necessary for a robust determi- , e . n m t, . ,. . ,. t ., . ■• 

3 hereafter Paper II). The most salient results of these studies can 

nation of the isochrone ages, which are very sensitive in cool , • ■ e ,, , s .. u . . ■ 

& ' 3 be summarized as follows: (a) active binaries are not as iron- 
stars to the choice of T e s, [Fe/H] and [a7Fe] . The lithium con- , c ■ . , , • ■ /UN ., . , 

L ' J L ' J dencient as previously claimed, (b) they present a large over- 

~Z ] ~Z '. rx, ., , ., abundance of most elements with respect to iron and (c) they 

Send offprint requests to: T. Morel, e-mail: r v ' J 

morel@astropa . unipa . it exhibit photometric anomalies leading to spuriously low (V-R) 

* Based on observations collected at the European Southern and ( V ~ ^ c °l° ur temperatures. 
Observatory, Chile (Proposals 64.L-0249 and 071. D-0260). These inferences were, however, somewhat tentative be- 

** Table |A~T1 is only available in electronic form at cause of the small number of objects studied. We have therefore 

http : //www . edpsciences . org/aa complemented these data by obtaining spectra for 8 additional 
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single-lined RS CVn binaries. Here we use the whole dataset 
to draw more robust conclusions regarding the behaviour of the 
photospheric abundances as a function of the activity level and 
to address the reliability of LTE abundance analyses at the up- 
per end of the stellar activity distribution. 

2. Observational material 

The spectra were obtained in July 2003 at the ESO 2.2-m 
telescope (La Silla, Chile) with the fiber-fed, cross-dispersed 
echelle spectrograph FEROS in the object+sky configuration. 
The spectral range covered is 3600-9200 A (resolving power 
/?~48 000), hence enabling a simultaneous coverage of Can 
H+K and the spectral lines used as abundance indicators. 

The 8 main programme stars (all of spectral type G8-K2 
IV— III) are drawn from the catalogue of chromospherically ac- 
tive stars of Strassmeier et al. (1993). Since fast rotators are 
not amenable to a curve-of-growth abundance analysis because 
of blending problems, only stars satisfying vsin/ ^ 10 km 
s were considered. Seven presumably single stars of similar 
spectral type, but with a much lower level of X-ray emission 
(Hiinsch, Schmitt & Voges 1998) were also observed. They 
were selected on the basis of their colours and absolute magni- 
tudes. Because of observational constraints, the reader should 
bear in mind that the control stars do not exactly match the 
properties of the target sample: they are on average slightly 
bluer and younger. Some basic properties of the targets and 
further details on the observing campaign (mean heliocentric 
Julian date of the observations, typical signal-to-noise ratio) 
can be found in Tables la-b. As will be apparent in the fol- 
lowing, one star originally classified as an active binary (HD 
28, 33 Psc) actually displays an extremely low level of activity 
(see also Walter 1985). 

Tasks implemented in the IRAF 1 package were used to 
carry out the standard reduction procedures for echelle spectra 
(i.e. bias subtraction, flat-field correction, removal of scattered 
light, order extraction and wavelength calibration). The con- 
tinuum normalization was performed using Kurucz synthetic 
spectra (Kurucz 1993) calculated for an initial guess of the at- 
mospheric parameters and metallicities. When available, tem- 
perature and metallicity values were taken from the compila- 
tions of Taylor (2003) and Cayrel de Strobel, Soubiran & Ralite 
(2001), respectively. These preliminary temperature estimates 
differ by less than 125 K from our final values derived from the 
excitation equilibrium of the iron lines. Otherwise, the temper- 
atures were determined from (B - V) indices (Alonso, Arribas 
& Martinez-Roger 1999, 2001) or, in a few cases, calibrated 
on the spectral type (Flower 1996). If unknown, solar metal- 
licity was assumed. We adopted for all the stars logg=2.7 cm 
s~ 2 . In the vast majority of cases, this estimate is in reson- 
able agreement both with previous works (Cayrel de Strobel 
et al. 2001) and with our final values. The vsin; values used 
are quoted in Tables la-b. The line-free regions were selected 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



from these synthetic spectra and then fitted by low-degree poly- 
nomials in the observed spectra. Every spectral line in our list 
was inspected in the 3 consecutive exposures obtained in order 
to identify the spectra with an imperfect continuum normal- 
ization or cosmic ray events. Any individual exposure affected 
by such problems was excluded during the co-adding proce- 
dure. Finally, we corrected for the star's systemic velocity by 
Doppler-shifting the spectra by the mean radial velocity of the 
lines under study. 

3. Analysis 

To ensure consistency, the abundance analysis carried out here 
is identical in all respects to the procedure employed in Paper 
II (to which the reader is referred to for full details). In par- 
ticular, we used the same line list (optimized for early K-type 
stars) and oscillator strengths. To calibrate the atomic data, we 
made use of a very high S/N FEROS moonlight spectrum and a 
Kurucz solar model with 7^=5777 K, log g=4.44 cm s~ 2 and 
a microturbulent velocity £=1.0 km s (Cox 2000). The oscil- 
lator strengths were adjusted until Kurucz's solar abundances 
(quoted in Table |A~T1 only available in electronic form) were 
reproduced. An extensive comparison between our log gf val- 
ues and literature data is presented in the Appendix of Paper II. 

The effective temperature and surface gravity were derived 
from the excitation and ionization equilibrium of the iron lines, 
while the microturbulent velocity was determined by requir- 
ing the Fe i abundances to be independent of the line strength. 
The abundances were derived from a differential LTE analysis 
using plane-parallel, line-blanketed Kurucz atmospheric mod- 
els (with a length of the convective cell over the pressure scale 
height, a-l/Hp-0.5, and no overshoot) and the current version 
of the MOOG software (Sneden 1973). The equivalent widths 
(EWs) were measured from Gaussian fitting and are quoted in 
Table lA. ll The lithium abundance was determined from a spec- 
tral synthesis of the Li i /16708 doublet (see Paper II). 

The reliability of abundance analyses relying on excitation 
temperatures may be questioned because of departures from 
LTE preferentially affecting the low-excitation iron lines (e.g. 
Ruland et al. 1980). Similarly to what has been done in Paper II, 
we have therefore repeated the whole procedure after rejection 
of the Fe i lines with a first excitation potential, ;f, below 3.5 eV. 
In this case, we used effective temperatures derived from T e g- 
(B - V) empirical calibrations (Alonso et al. 1999). As will be 
shown below, these values turn out in most cases to be slightly 
lower than the temperatures derived from the excitation equi- 
librium of the Fe i lines. The resulting abundance ratios will 
not be considered further, as they are identical, within the un- 
certainties, to the values derived using our main method. 

4. Results 

The atmospheric parameters (Teg, logg and £) and abundances 
are given in Tables la-b and 2a-b, respectively. Note that 
the abundances of Li, Na, Mg and Ca have been converted 
into NLTE values. The corrections range from Ae=log(e)NLTE- 
log(e) LTE =+0. 10 to +0.26 for Li i ,16708 (Carlsson et al. 1994), 
from -0.09 to +0.06 dex for Nai /16154 and from +0.04 to 
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Fig. 1. Differences between the excitation and colour temperatures, as a function of the activity indices R' nK (top panels), Rx 
(middle panels) and maximum amplitude of the wave-like photometric variations in V band, AV wave (bottom panels). The varia- 
tions are shown for temperatures derived from (B - V) (left-hand panels), (V - R) (middle panels) and (V - I) data (right-hand 
panels). The active binaries and stars in the control sample are indicated by filled and open circles, respectively (in the online 
version of this journal, red and blue circles refer to active and control stars, respectively). The crosses are upper limits. 



+0.15 dex for Mgi ,15711 (Gratton et al. 1999). For conve- 
nience, these corrections are listed on a star-to-star basis in 
Table IB. II (we also include the stars analyzed in Paper II). 
For the Ca lines, we used for all the stars Ae=+0.09, +0.07 
and +0.01 for Cai ,16166, 6456 and 6500, respectively (Drake 
1991). For the other chemical species, no NLTE corrections for 
the relevant range of physical parameters are available in the 
literature. 

The radiative loss in the Can H+K lines in units of the 
bolometric luminosity, R' nK , was used as primary activity in- 
dicator. This quantity was derived from our spectra following 
Linsky et al. (1979). The relative emission-line fluxes of the 
Ca ii H+K lines were calculated by integrating the unnormal- 
ized spectra in the wavelength range encompassing the emis- 
sion components, and subsequently dividing these quantitities 
by the total counts integrated in the 3925-3975 A bandpass. 
The conversion to absolute fluxes was performed using an em- 
pirical calibration between (V - R) and the absolute surface 
flux in the 3925-3975 A wavelength range, T^o (Linsky et al. 
1979). Finally, these two quantities were corrected for photo- 
spheric contribution, and the sum divided by the bolometric 



luminosity to obtain R' HK (see Paper II for further details). The 
colour anomalies discussed in Sect. 15.11 lead to a systematic, 
dramatic underestimation of ^50, and ultimately R' nK (up to 0.5 
dex), when adopting the observed (V - R) colour. We there- 
fore used instead the colour appropriate for a star with the de- 
rived T eff and [Fe/H] (eq. (5) of Alonso et al. 1999). We re- 
gard this approach as a major improvement compared to what 
has been done in Paper II. Consequently, the activity indices 
for the 6 stars discussed in this paper have been re-calculated 
accordingly. The bolometric luminosities were estimated from 
theoretical isochrones (see Sect l5.2t . We also define Rx, which 
is given as the ratio between the X-ray (in the ROSAT [0.1- 
2.4 keV] energy band) and bolometric luminosities. Finally, we 
use the maximum amplitude of the wave-like photometric vari- 
ations in V band, AV wave , as a proxy of stellar spottedness (val- 
ues from Strassmeier et al. 1993). The activity indices, as well 
as the source of the X-ray data, are given in Tables la-b. 

The kinematic data for the programme stars were calcu- 
lated following Paper II (Tables la-b). The control stars have 
kinematical properties typical of the thin disk population, but 
some RS CVn binaries have large rotational lags (V~-50 km 
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Fig. 2. Positions of the programme stars in HR diagrammes: (a) HD 218527, (b) HD 4482, (c) HD 182776, (d) HD 217188, (e) 
HD 28, if) HD 181809, (g) HD 204128, (h) HD 154619, (/) HD 1227, (/') HD 205249, (*) HD 211391, (0 HD 156266, (m) HD 
17006, («) HD 19754 and (o) HD 202134. The metallicity of the isochrones refers to a solar iron abundance: log e o (Fe)=7.50 (Yi 
et al. 2001 ; Kim et al. 2002). The age of the isochrones (in Gyr) is indicated in the left-hand panels. Symbols as in Fig.^ 



s _1 ) pointing to a possible thick disk membership. We follow 5. Reliability of the atmospheric parameters 



Bensby, Feltzing & Lundstrom (2004) to estimate the proba- 
bility that a star belongs to the thin/thick disk or to the halo, 
using the calculated velocity components with respect to the 
LSR (U, V and W) and considering the characteristic velocity 
dispersions of these stellar populations (cry, cry and C7>). It is 
assumed that the galactic space velocities in the solar neigh- 
bourhood can be approximated by a Gaussian distribution. The 
other quantities entering in the calculations are the observed 
local fraction of stars, X, and the asymmetric drift, y asym , cor- 
responding to each population. All the relevant values are taken 
from table 1 of Bensby et al. (2004). The likelihood for a thin or 
thick disk membership can be parameterized by the ratio of the 
respective probabilities for both components, D/TD (Bensby 
et al. 2004). Only HD 83442 is very likely to be a thick disk 
star (D/TD-0.02, i.e. the star is 50 times more likely to belong 
to the thick disk than to the thin disk), whilst the case of HD 
10909 and HD 119285 is more ambiguous (D/TD=0.S6 and 
0.60, respectively). These 3 active binary systems have been 
discussed in Paper II. All the other stars in our sample are very 
likely to belong to the thin disk (D/TD > 8). None is part of the 
galactic halo. 



An abundance analysis of active stars relying on effective tem- 
peratures and surface gravities derived from classical 1-D LTE 
atmospheric models might be questioned because of the exis- 
tence of large-scale photospheric structures (e.g. plages, cool 
spots) and an overlying, prominent chromosphere. Before dis- 
cussing the abundance patterns, let us first examine the relia- 
bility of our derived atmospheric parameters. In the following 
discussion, we shall combine our results with data taken from 
Paper II for 6 additional RS CVn binaries. 

5.1. Effective temperatures 

A consistency check on our temperature estimates is provided 
by a comparison with colour temperatures computed from the 
(B - V), (V - R) and (V - /) indices, using calibrations based 
on the infrared flux method (Alonso et al. 1999). The photo- 
metric data used, as well as the literature sources are quoted 
in Tables la-b. Preference was given to repeated observations 
in order to minimize possible colour variations arising from 
the rotational modulation of active regions. In most cases, the 
adopted colours are the mean of more than 15 (and up to 70) 
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Table la. Physical parameters of the active binaries. Spectral type, rotational and orbital periods from Strassmeier et al. (1993). 
We quote the typical S/N ratio at 6700 A after combination of the individual exposures. The A{V) values are derived from the 
empirical galactic extinction model of Arenou, Grenon & Gomez (1992). The colours in the Cousins and Kron photometric 
systems have been transformed into Johnson values following Bessell (1979). AV wave is the maximum amplitude of the wave-like 
photometric variations in V band (Strassmeier et al. 1993). v r is the radial velocity adopted in the computation of the velocity 



components. The peculiar space velocity, S , is given by: 
masses derived from the evolutionary tracks. 



S = (U + V + W ) ' . The last two rows give the stellar ages and 
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time-resolved measurements. In case of a metallicity depen- 
dent relation, we used our derived [Fe/H] values. The excitation 
temperatures of the active binaries appear marginally higher 
than the values derived from T e ff-(B - V) calibrations: <r exc - 
7 , C oiour>=+80±46 K. This fair agreement is gratifying, but we 
cannot rule out that it is circumstantial in view of previous 
claims of a (B - V) excess in active stars (e.g. Amado 2003; 
Stauffer et al. 2003). We also warn the reader that our excita- 
tion temperatures may be affected by other shortcomings (e.g. 
NLTE effects). In sharp contrast, the mean difference amounts 
to +201+36 and +215+35 K for (V-R) and (V-I), respectively. 
As discussed in Papers I and II, the differences are unlikely to 
be accounted for by spottedness or by the existence of a faint, 
cooler companion. 

We show in Fig.^a comparison between the excitation and 
colour temperatures, as a function of the three activity indices 
defined above. There is a general trend for an increasing dis- 
agreement with higher R' HK and Rx values. The false alarm 
probabilities for a chance correlation between the temperature 
differences and the activity indices are quoted in Table[5] An in- 
creasing discrepancy between the excitation temperatures and 



the (V - R)- and (V - 7)-based estimates at higher activity level 
is quantitatively confirmed for the Rx and, to a lesser extent, 
for the R' nK data. This conclusion holds irrespective of whether 
the whole dataset or only the active binaries are considered. 
There is a hint of a similar trend for the (B - V) data, but the 
correlation is not statistically significant. 

The tendency for the (V - R) and (V - I) indices to yield 
systematically lower temperatures than {B - V) has been re- 
ported in several occasions for active stars, either in the field 
(Fekel, Moffett & Henry 1986) or in young open clusters 
(Garcia Lopez, Rebolo & Martin 1994; Randich et al. 1997; 
Soderblom et al. 1999). Although early attempts to link these 
colour anomalies to activity have been inconclusive in disk 
stars, presumably because of a limited range in activity level 
(Rucirisky 1987), evidence for a (V-I) excess in the most active 
Pleiades stars has been presented by King, Krishnamurthi & 
Pinsonneault (2000). A distortion of the spectral energy distri- 
bution because of activity-related processes is indeed expected, 
but the quantitative effect on the broad band colours is diffi- 
cult to assess (e.g. Stuik, Bruls & Rutten 1997). Chromospheric 
models for M-type dwarfs indicate a substantial modulation, as 
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Table lb. Physical parameters of the control stars. 





HD 1227 


HD 4482 


HD 17006 


HD 154619 


HD 156266 


HD 21 1391 


HD 218527 


Spectral Type 


G8 II-III 


G8II 


Kim 


G8 III-IV 


K2III 


G8 III-IV 


G8 III-IV 


vsini (km s -1 ) 


<i.oi') 


si.o") 




1.3W 


<15t!2] 


<15l' 2 ) 


3.0 [I3] 


5 IN 


190 


205 


190 


190 


180 


225 


115 


A(V) (mag) 


0.070 


0.069 


0.052 


0.247 


0.267 


0.096 


0.070 


(B - V) (mag) 


0.899 1141 


0.949115] 


0.864 1141 


0.795 [16] 


1.069 [14] 


0.96 1 [14] 


0.889 [141 


(V-/?)o(mag) 


0.683 [17] 


0.733 [17] 


0.677 [I7] 




0.734 [6] 


0.676 [6] 




(V-/)o (mag) 


U44 [n] 


1.204 [17] 


1.113 [17] 




1.251 [6] 


1.130 [6] 




rcoiou, (B - V) (K)" 


5076 


4922 


5224 


5295 


4739 


4955 


5006 


T colom (V-R)(K) b 


4997 


4837 


5034 




4853 


5023 




T aiom (V - I) (K) fc 


4985 


4872 


5046 




4788 


5012 




Texc (K) 


5115+67 


4995±67 


5332+37 


5260+44 


4665+58 


5025+60 


5040+88 


log g (cm s~ 2 ) 


3.04±0.10 


2.88±0.08 


3.80+0.09 


3.38+0.09 


2.50+0.13 


2.84+0.08 


2.79+0.20 


f (kms- 1 ) 


1.33±0.06 


1.47±0.07 


1.22+0.06 


1.28+0.07 


1.48+0.06 


1.47+0.07 


1.23+0.12 


log (R' HK f 




-4.23±0.10 


^t.31+0.10 


^t.52+0.11 


-5.05+0.12 






log(L x ) (ergs s l ) d 


29.39±0.23 [18] 


29.83+0. 10 [18] 


29.45±0.07 [18] 


29.52±0.20 [18] 


28.97+0.18 1181 


29.42+0. 13 1181 


29.86+0. 14 [181 


log (RxY 


-5.92+0.23 


-5.34±0.10 


^t.83+0.07 


-5.72+0.20 


-6.46+0.18 


-6.00+0.13 


-5.39+0.14 


v r (km s" 1 ) 


-0.3±0.2 [11 


-5.1 + 1.51') 


+ 13.0+2.0 1151 


-22.8±0.2 [1] 


-2.3±0.9 [19] 


-14.7+0.9 1191 


-17.8±2.0 [19] 


[/(kins" 1 ) 


8.05±2.56 


-21.0+1.0 


2.63+0.50 


-27.9+0.4 


^1.06+0.96 


-40.8+1.2 


-83.6+7.8 


Vfkms 4 ) 


-3.13+0.77 


-27.7+1.6 


-23.3+1.2 


-0.90+1.61 


-27.1 + 1.2 


-24.2+0.9 


-1.96+2.19 


W (kins' 1 ) 


-9.95±0.59 


-12.0+1.3 


-13.0+1.8 


-33.0+1.6 


-8.83+0.56 


-16.7+1.2 


6.80+1.76 


5 (kms- 1 ) 


13.2+1.6 


36.8+1.4 


26.8+1.3 


43.2+1.3 


28.8+1.2 


50.3+1.1 


83.9+7.8 


T is0 (Gyr) 


5 +01 
"• J -o.i 


0-811 


9 q+o.l 

^•'-o.i 


7 +al 
"•'-o.i 


6 +01 
"•"-o.i 


40 +005 

"•^"-0.05 


6 +0 2 
"•"-o.i 


M(Mq) 


-2.8 


-2.4 


-1.4 


-2.6 


-2.6 


-3.0 


-2.6 



Key to references — [1] de Medeiros & Mayor (1999); [2] Randich, Gratton & Pallavicini (1993); [3] Walter (1985); [4] Evans & Koen 

(1987); [5] Hipparcos data; [6] Johnson et al. (1966a); [7] Dempsey et al. (1993); [8] Schwope et al. (2000); [9] Neuhauser et al. (2000); [10] 

Dempsey et al. (1997); [11] Balona (1987); [12] Herbig & Spalding (1955); [13] Gray & Nagar (1985); [14] Johnson et al. (1966b); [15] 

Haggkvist & Oja (1987); [16] Oja (1991); [17] Eggen (1989); [18] Hiinsch et al. (1998); [19] SIMBAD database. 

The typical 1-cr uncertainty is 150 K. 

b The typical 1-cr uncertainty is 100 K. 

c A blank indicates that the Can H+K profiles were not filled in by emission. The errors take into account the uncertainties in (V - R)o and T ef j. 

d All luminosities were rescaled to the Hipparcos distances. The errors take into account the uncertainties in the X-ray counts and in the 

distances. 



a function of the activity level, not only of the emergent con- 
tinuum flux in the UV domain, but also in the optical and near- 
infrared wavebands (Houdebine et al. 1996; see, however, West 
et al. 2004). On the observational side, the infrared excess in- 
trinsic to some RS CVn binaries has been claimed to diagnose 
the existence of cool, circumstellar material (e.g. Scaltriti et al. 
1993, and references therein). 

5.2. Surface gravities 

We show in Fig.|2]the position of our stars in HR diagrammes, 
along with theoretical isochrones chosen to match as closely 
as possible our derived [Fe/H] and [a/Fe] values (Yi et al. 
2001; Kim et al. 2002). We use Hipparcos data and the exci- 
tation temperatures. The stellar ages and masses derived from 
the evolutionary tracks are given in Tables la-b. The surface 
gravities determined from the ionization equilibrium of the Fe 
lines appear systematically lower for the RS CVn binaries than 
the values derived from the isochrones (<loggj on i-log,gi so >=- 
0.21+0.06 dex). In contrast, there is a good agreement for the 
control stars: <loggi on i-log,gi so >=-0.01+0.08 dex. This dis- 
crepancy in the gravity estimates could be plausibly caused 
by an overionization of iron (see below). There is a hint 



that the surface gravities determined from fitting the wings of 
collisionally-broadened lines tend to be systematically higher 
than values derived from the ionization equilibrium of the Fe 
lines (Paper I). Other interpretations are, however, possible. 
Because T e g and logg are interdependent, for instance, such 
a gravity offset could be an artefact of a small, systematic un- 
derestimation of the effective temperatures (at the ~70 K level). 
Summarizing, we conclude that our T e s and logg esti- 
mates are liable to systematic errors of the order of 150 K and 
0.25 dex, respectively. The sensitivity of the abundance ratios 
against such changes is examined in Tabled The impact is lim- 
ited and is typically comparable to the abundance uncertainties 
(A[M/Fe] < 0.12 dex). 



6. Abundance patterns 

As can be seen in Fig. [5] there is a fairly clear dichotomy be- 
tween the abundance patterns of the two samples of active and 
inactive stars. The active binaries exhibit as a class a chemical 
composition departing more conspicuously from the solar mix 
(compare the left- and right-hand panels). 

Figure 0] shows the abundance ratios as a function of 
[Fe/H]. The increasing abundance ratios of several chemical 
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Fig. 3. Abundance patterns of the RS CVn binaries (left-hand panels) and the stars in the control sample (right-hand panels). The 
chemical elements are identified in the right-hand, upper panel. The position of barium has been shifted for the sake of clarity to 

Z=33. 



species (e.g. the a elements) with decreasing [Fe/H] is a well- 
known feature of galactic disk stars, and may be naturally 
explained in terms of a varying relative contribution of Type 
la/Type II supernova ejecta to the chemical composition of the 
natal interstellar cloud. However, the large overabundances ob- 
served cast some doubt on this interpretation. We overplot in 
Fig. |3 the characteristic abundance patterns of kinematically- 
selected samples of thin and thick disk FG dwarfs (Bensby et 
al. 2003). The overabundances presented by the active binaries 
appear generally much more marked at a given metallicity than 
for either of these two populations. For instance, [Ca/Fe]~0.35 
dex at [Fe/H]=-0.4, while [Ca/Fe]~0.1 and 0.2 dex only for 
thin and thick disk stars, respectively (we recall that very few 
stars in our sample are expected to belong to the thick disk; 
Sect. |3. The dramatic Al overabundances at near solar metal- 
licities (up to 0.8 dex) are also completely at variance with what 



is expected for stars in the solar neighbourhood. The origin of 
the high Na and Ba abundances in some control stars is unclear, 
but might be related to systematic errors in the atomic data. The 
abundances are derived from a single line in both cases. 



The RS CVn systems often display higher abundance ratios 
than the control stars at a given [Fe/H] value (Fig.@J- To test the 
hypothesis that it is related to differences in the activity level, 
we plot in Fig.[5]the abundance ratios as a function of R' nK and 
^x- The positive correlation seen for several elements (e.g. Al, 
Ca) suggests that the observed overabundances may arise from 
activity-related processes. This hypothesis is examined in the 
following. 
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Table 2a. Abundances of the active binaries. We use the usual notation: [A/B]=log [Af(A)/Af(B)] + -log [N(A)/N(B)]q. The 
corresponding l-cr uncertainties and the number of lines used, N, are indicated. Blanks indicate that no EWs for the element 
in question could be reliably measured. The last row gives the lithium abundances. In this case we use the usual notation: 
log e(Li) = log A/(Li) - 12.0. The abundances of Na, Mg, Ca and Li have been corrected for departures from LTE (see text and 
TablelinV 





HD28 


HD 19754 


HD 181809 


HD 182776 


HD 202134 


HD 204128 


HD 205249 


HD 217188 


Name 


33Psc 


EL Eri 


V4138 Sgr 


V4139 Sgr 


BNMic 


BHInd 


AS Cap 


AZPsc 


[Fe/H] 


0.01±0.08 


-0.39+0.08 


-0.09+0.04 


-0.06+0.16 


-0.06+0.07 


-0.02+0.11 


0.12+0.09 


-0.16+0.07 


N 


47 


37 


38 


20 


45 


30 


36 


43 


[Na/Fe] 


0.13+0.05 


0.50+0.08 


0.30+0.04 


0.73+0.25 


0.53+0.04 


0.52+0.13 


0.56+0.05 


0.37+0.04 


N 


1 


1 


1 


1 


1 


1 


1 


1 


[Mg/Fe] 


0.08±0.05 


0.26+0.05 


0.15+0.03 




0.16+0.04 


0.08+0.07 


0.09+0.04 


0.08+0.03 


N 


1 


1 


1 




1 


1 


1 


1 


[Al/Fe] 


0.13+0.05 


0.35+0.04 


0.21+0.04 


0.82+0.11 


0.37+0.05 


0.36+0.05 


0.30+0.03 


0.27+0.03 


N 


2 


1 


2 


1 


2 


1 


2 


2 


[Si/Fe] 


0.12+0.08 


0.03+0.08 


-0.01+0.06 


0.05+0.14 


0.07+0.10 


-0.04+0.09 


0.08+0.10 


-0.03+0.06 


N 


9 


5 


6 


4 


6 


4 


7 


5 


[Ca/Fe] 


0.11 ±0.07 


0.32+0.06 


0.22+0.04 


0.14+0.09 


0.27+0.05 


0.25+0.14 


0.18+0.05 


0.22+0.04 


N 


3 


3 


3 


3 


3 


3 


2 


3 


[Sc/Fe] 

N 

[Ti/Fe] 


0.11+0.06 

1 

0.02+0.04 
















0.17+0.05 


0.14+0.03 




0.18+0.04 


0.15+0.09 


0.26+0.09 




N 


1 


1 


1 




1 


1 


1 




[Cr/Fe] 


-0.10+0.10 


0.10+0.11 


0.11+0.05 


0.02+0.21 


0.14+0.05 


0.14+0.12 


0.17+0.16 


0.03+0.09 


N 


3 


2 


1 


2 


1 


2 


3 


3 


[Co/Fe] 




0.05+0.08 


0.03+0.04 




0.10+0.05 


0.17+0.13 


0.19+0.05 




N 




1 


1 




1 


1 


1 




[Ni/Fe] 


0.01+0.07 


-0.10+0.07 


-0.07+0.05 


-0.12+0.17 


-0.03+0.08 


-0.03+0.06 


0.01+0.09 


-0.13+0.06 


N 


10 


8 


10 


5 


10 


8 


10 


8 


[Ba/Fe] 


0.10+0.08 


0.36+0.09 


0.29+0.06 


-0.15+0.27 


0.15+0.06 


0.12+0.15 


0.06+0.06 


0.27+0.08 


N 


1 


1 


1 


1 


1 


1 


1 


1 


loge(Li) 


<0.87 


1.20+0.13 


1.23+0.09 


1.18+0.22 


0.95+0.15 


1.49+0.17 


2.03+0.10 


1.22+0.12 



7. Discussion 



Cool spot groups are known to cover a large fraction of the 
stellar photosphere in chromospherically active stars. For one 
star in our sample (HD 1 19285; see Paper II), Saar, Nordstrom 
& Andersen (1990) derived a spot filling factor in the range 
0.18-0.37 from a modelling of the TiO bandhead at 8860 A. 
To investigate the impact of these features on our results, we 
have created a composite Kurucz synthetic spectrum of three 
fictitious stars with r e ff=4830 K, covered by cooler regions 
(Ar e ff=-1000 K) with different covering factors, / s =0, 30 and 
50% (Paper II, see also Neff, O'Neal & Saar 1995). We then 
derived the abundances from these composite spectra follow- 
ing exactly the procedure used for the analysis of the active 
binaries. The differences between the results obtained for the 
unspotted star are given in Table [5] As can be seen, the exis- 
tence of cool spots generally leads to a systematic overestima- 
tion of the abundance ratios with respect to iron. The impact on 
the resulting abundances can be significant and may account 
for the overabundances observed for some elements (e.g. Ca). 
This mostly results from a filling in of the iron line profiles, 
a behaviour which seems supported by the apparent iron de- 
ficiency at high activity level seen in Fig. [6] (see also Cayrel, 



Cayrel de Strobel & Campbell 1985; Drake & Smith 1993). 2 
If this interpretation is correct, an asymmetrical distribution of 
cooler photospheric regions is expected to lead to a variation 
of the derived abundances as a function of the rotational phase. 
Ottmann, Pfeiffer & Gehren (1998) carried out time-resolved 
spectroscopy of several RS CVn binaries, but did not find con- 
vincing evidence for such an effect. 

Most spectral lines used in this study are generally strong 
and of low excitation, and are therefore bound to form in upper 
photospheric layers where the effect of the low chromosphere 
may be appreciable. We have investigated in Paper II the in- 
fluence of a chromospheric temperature rise on the abundance 
ratios. To this end, a LTE abundance analysis was carried out 
using modified Kurucz atmospheric models characterized by a 
temperature gradient reversal in the outermost regions and an 
overall heating (up to 190 K) of the deeper photospheric lay- 



2 We note, however, that the two samples of active and inactive stars 
are likely drawn from two distinct populations in terms of stellar ages. 
The median evolutionary age of the control sample is only about 0.6 
Gyr, against 2.4 Gyr for the active binaries. Based on the statistical 
relationship in galactic disk stars between the stellar age and the iron 
content (Rocha-Pinto et al. 2000), we might expect a ~0.2 dex offset 
in the mean metallicity of the two samples. 
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HD 1227 


HD 4482 


HD 17006 


HD 154619 


HD 156266 


HD 211391 


HD 218527 


[Fe/H] 


0.19+0.07 


0.05+0.07 


0.36±0.06 


0.13+0.06 


0.25±0.08 


0.23±0.07 


-0.08±0.09 


N 


40 


37 


51 


48 


35 


39 


36 


[Na/Fe] 


0.18±0.04 


0.30±0.05 


0.12+0.05 


0.02±0.04 


0.30±0.05 


0.28±0.04 


0.11+0.14 


N 


1 


1 


1 


1 


1 


1 


1 


[Mg/Fe] 


-0.05±0.04 


0.04±0.05 


-0.04±0.03 


-0.04±0.03 


0.01±0.04 


-0.05+0.04 


-0.09±0.07 


N 


1 


1 


1 


1 


1 


1 


1 


[Al/Fe] 


0.00±0.03 


0.13±0.03 


0.07±0.06 


0.02+0.05 


0.24±0.05 


0.04±0.04 


0.02+0.08 


N 


2 


2 


2 


2 


2 


2 


1 


[Si/Fe] 


0.04±0.06 


0.08±0.08 


0.05+0.07 


~0.03±0.04 


0.16±0.04 


0.08±0.07 


0.04±0.13 


N 


6 


7 


7 


6 


3 


7 


7 


[Ca/Fe] 


0.09±0.04 


0.13±0.07 


0.02±0.05 


0.10±0.05 


0.09±0.10 


0.09±0.05 


0.06±0.06 


N 


3 


3 


3 


3 


3 


3 


3 


[Sc/Fe] 


-0.03±0.07 






0.03±0.06 


0.05±0.07 






N 


1 






1 


1 






[Ti/Fe] 


-0.04±0.05 


0.07+0.04 


0.02+0.07 


0.10±0.09 


0.05+0.04 


0.00±0.05 


0.08±0.13 


N 


1 


1 


1 


1 


1 


1 


1 


[Cr/Fe] 


-0.08±0.05 


-0.03±0.05 


0.02+0.08 


-0.01 ±0.05 


0.05±0.09 


0.00±0.05 


-0.04±0.15 


N 


1 


1 


2 


1 


2 


1 


1 


[Co/Fe] 


0.02±0.05 




0.14±0.05 


0.00±0.05 








N 


1 




1 


1 








[Ni/Fe] 


-0.03±0.06 


-0.06±0.08 


0.05+0.06 


-0.09±0.05 


0.05±0.09 


-0.01±0.07 


-0.06±0.08 


N 


10 


10 


10 


10 


10 


10 


8 


[Ba/Fe] 


0.30±0.06 


0.25+0.06 


-0.09±0.07 


0.40±0.05 


-0.07±0.08 


0.23+0.05 


0.25+0.13 


N 


1 


1 


1 


1 


1 


1 


1 


loge(Li) 


1.16+0.18 


1.31±0.12 


<1.18 


<1.03 


<0.75 


<0.95 


<0.51 



Table 3. False alarm probabilities for a chance correlation between the temperature differences and the activity indices (all values 
expressed in per cent). The numbers in brackets refer to the active binary dataset alone. 





T^-T(B - V) 


T CXC -T(V-R) 


T^-T(V-I) 


log^W 


13.7 (58.2) 


2.9 (3.3) 


5.3 (5.8) 


log (Rx) (whole dataset)" 


7.4(95.1) 


0.2 (4.6) 


0.1 (1.2) 


log (R x ) (detections only) 


12.7(95.1) 


0.3 (4.6) 


0.3 (1.2) 


^ 'wave 


(41.0) 


(19.9) 


(8.9) 



" We used the Kendall's t method in the case of censored data (HD 28 is undetected in X-rays). 



ers (see fig. 3 of Paper II). We considered two model chromo- 
spheres meant to be representative of K-type giants (Kelch et 
al. 1978) and RS CVn binaries (Lanzafame, Busa & Rodono 
2000). Excluding Sc and Ba, the incorporation of a chromo- 
spheric component yields abundance ratios that do not differ 
by more than 0.13 dex from those given by classical Kurucz 
models. Taking into account chromospheric heating may have 
a significant impact on the derived abundances of some ele- 
ments, but appears insufficient to account for the anomalously 
high Na and Al abundances, for instance. 

Another important, related issue is the fact that the NLTE 
corrections applied to the abundances of some elements (Na, 
Mg, Ca) are based on 'conventional' atmospheric models not 
including a chromospheric component, the latter being a key 
ingredient in NLTE line formations studies (e.g. Bruls, Rutten 
& Shchukina 1992). In a companion paper, we have recently 
discussed the behaviour of the optical oxygen abundance in- 
dicators in active stars (Morel & Micela 2004). While [Oi] 
/16300 appears to be largely unaffected by chromospheric ac- 



tivity, the 7774 A-0 1 triplet yields abundances dramatically in- 
creasing with R' nK or Rx, even after 'canonical' NLTE correc- 
tions are applied. These two diagnostics have about the same 
sensitivity (in an absolute sense) against atmospheric parame- 
ter changes, but a completely distinct behaviour in terms of de- 
partures from LTE (the O i triplet is strongly affected, whereas 
[O i] /(6300 is insensitive to these effects). This discrepancy be- 
tween the two indicators may therefore indicate that the NLTE 
corrections to apply to the permitted oxygen lines in chromo- 
spheric ally active stars largely exceed what is anticipated for 
their inactive analogues. It is conceivable that a similar con- 
clusion, albeit to a lesser extent, also holds for other chemical 
species. As far as iron is concerned, however, the fair agree- 
ment between the excitation/(B-V)-based temperatures or the 
ionization/isochrone gravities (Sect.|5j does not argue in favour 
of strong deviations from LTE in our sample. One may ex- 
pect overionization to be an important process in active stars 
in view of their strong UV excess (Stauffer et al. 2003), and 
to manifest itself by an overall depletion of the neutral ioniza- 
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Fig. 4. Abundance ratios as a function of [Fe/H]. The active binaries and stars in the control sample are indicated by filled 
and open circles, respectively. The possible thick disk stars (HD 10909, HD 83442 and HD 119285; see text) are indicated by 
filled triangles. Colour coding as in Fig.^ We define the mean abundance ratio of the a-synthezised elements, [a/Fe], as the 
unweighted mean of the Mg, Si, Ca and Ti abundances. The thick dashed and dotted lines show the characteristic trends of 
kinematically-selected samples of thin and thick disk stars, respectively (Bensby et al. 2003), The data for Sc, Co and Ba are 
taken from Reddy et al. (2003). 



tion stage. Evidence for an unexpectedly large overionization 
of iron in young open cluster stars is indeed recently emerg- 
ing (e.g. Schuler et al. 2003; Yong et al. 2004). If explained in 
terms of NLTE effects, the large overabundances derived from 
neutral lines (e.g. Nai, Ah) would apparently require a com- 
peting physical mechanism capable of dramatically overpopu- 
lating the neutral stage in such a way that it overwhelms the 
action of overionization (see, e.g. Bruls et al. 1992 for a dis- 
cussion of the various NLTE mechanisms possibly at work). 

Contrary to the situation for some globular cluster stars 
ascending the red giant branch, it is unlikely that the surface 
abundance anomalies observed result from deep mixing pro- 
cesses not accounted for by standard evolutionary models. One 



could speculate that this extra mixing mechanism is related in 
some way to the spinning up of the stars in asynchronous X- 
ray binaries because of the transfer of orbital angular momen- 
tum (see, e.g. Denissenkov & Tout 2000). This hypothesis can- 
not be discarded on the basis of the generally high Li abun- 
dances derived (Table l2al Paper II), as little Li depletion might 
not provide compelling evidence against prior mixing phases 
(Charbonnel & Balachandran 2000). 3 This extra mixing leads 
in some red giants to an enhancement of the O and Mg con- 



3 No relationship was found in our sample between the lithium con- 
tent and various quantities, including the evolutionary age, the activity 
level or the rotational velocity (see Paper II in the latter case). 
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Fig. 5. Abundance ratios as a function of the activity indices R' HK (upper panels) andRx (bottom panels). Symbols as in Fig.|4] 



tent accompanied by a decrease of the Al and Na abundances 
(e.g. Kraft et al. 1997), but these observational signatures are 
not seen in our sample. 

Independently from the various activity processes dis- 
cussed above, we call attention to the fact that temperature ef- 
fects may also play a significant role in inducing abundance 
peculiarities. In close analogy to the behaviour of the oxygen 
triplet (Morel & Micela 2004), the abundance ratios of several 
elements are antic orrelated with T e g (Fig. 0. Similar trends 



for some chemical species, but over a wider temperature range 
(4600-6400 K), have been reported in a sample of planet-host 
stars by Bodaghee et al. (2003). The dependence with respect 
to T e ff is qualitatively different, however, being even opposite 
in some cases to what we observe (e.g. Ca). Since the coolest 
stars in our sample (i.e. the RS CVn binaries) are also the most 
active, it is very difficult to disentangle temperature and ac- 
tivity effects. Unfortunately, few inactive stars with near-solar 
metallicity and in the relevant temperature range (T^ 5= 5000 
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Table 4. Effect of changes in the stellar parameters on the abun- 
dance ratios derived for HD 10909 (see Paper II). The iron con- 
tent is given by the Fe i lines. 
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K) have abundances reported in the literature (e.g. Feltzing & 
Gustafsson 1998; Affer et al. 2004; Allende Prieto et al. 2004 
for the most recent works). A systematic, homogeneous abun- 
dance study of a large sample of inactive K-type stars seems 
worthwhile for a proper interpretation of Fig.Q Of relevance is 
the substantial overionization of some metals recently reported 
in K-type field stars (e.g. Allende Prieto et al. 2004). This devi- 
ation from ionization equilibrium may indicate caveats in our 
modelling of the atmospheres of cool stars or, perhaps more 
likely, a still fragmentary understanding of line formation in 
these objects. In particular, the increasing magnitude of overi- 
onization effects with decreasing temperature is not predicted 
by current NLTE calculations (e.g. Thevenin & Idiart 1999). 

8. Conclusion 

The analysis of this enlarged sample of RS CVn binaries has 
allowed us to place the conclusions presented in Paper II on a 
much firmer footing. In particular, we confirm the existence in 
chromospherically active stars of a (V - R) and (V - I) excess 
leading to spuriously low colour temperatures. 

The active binary systems analyzed are characterized by 
abundance ratios generally at odds with the global patterns ex- 
hibited by galactic disk stars. We have argued that systematic 
errors in the determination of the atmospheric parameters are 
not expected to have a very significant impact on the abundance 
patterns (Table|4}. The effect of cool spot groups is to lead to a 
systematic overestimation of the abundance ratios with respect 
to iron (Table[3}, in qualitative agreement with the general pat- 
tern observed. As such, this constitutes an attractive explana- 
tion for the abundance peculiarities discussed in this paper. The 
observed correlation between the abundance ratios of some el- 
ements and the stellar activity level (Fig. |5} may also fit into 
this picture, considering the statistical relationship between the 
incidence of starspots and the X-ray emission (Messina et al. 
2003). However, the lack of an a priori knowledge of, e.g. 
the temperature difference between these cool features and the 



Table 5. Effect of cool spots on the derived stellar parameters 
and abundances. The quoted differences are given with respect 
to the results obtained for a Kurucz synthetic spectrum calcu- 
lated for the physical parameters and iron content of HD 10909 
(see Paper II). 
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undisturbed photosphere or the spot covering factor makes dif- 
ficult a quantitative foray into the relevance of this scenario. In 
addition, the effect of other competing mechanisms (e.g. plage 
areas, chromospheric heating) should, in principle, also be si- 
multaneously taken into account. The anomalous abundances 
presented by some elements (e.g. Na, Al) may betray the action 
of an additional process causally linked to chromospheric activ- 
ity. Unexpectedly large departure from LTE arising from non- 
thermal line excitation is a plausible candidate, but the physi- 
cal process capable of inducing such an overpopulation of the 
neutral stage remains to be identified. Overionization would be 
instead expected to play a significant role in stars with such a 
UV excess. 

Finally, it cannot be ruled out that general limitations in 
our current understanding of K-type stars with extended atmo- 
spheres (NLTE line formation, atmospheric modelling; Fields 
et al. 2003) are at least partly responsible of the odd abundances 
observed (Fig.0. 

Determinations of the coronal abundances in chromospher- 
ically active stars have flourished in recent years (e.g. Sanz- 
Forcada, Favata & Micela 2004) and have given an impetus 
to studies of chemical fractionisation processes in stellar coro- 
nae (e.g. Drake 2003). Unfortunately, the lack of photospheric 
abundance determinations for elements other than iron (which 
results in the customary use of the solar pattern as reference) 
may lead to misleading results and impedes a clear interpreta- 
tion of this phenomenon. Although our initial goal when this 
project was initiated was to fill this caveat by determining ac- 
curate photospheric abundances for chemical species of prime 
interest in the context of X-ray studies (namely, O, Al, Ca, Ni, 
Mg, Fe, Si), 4 we have shown here that the photospheric abun- 



4 Carbon and nitrogen spectral lines are not measurable in our spec- 
tra owing to their intrinsic weakness and severe blending. 
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Fig. 6. Iron content as a function of the activity indices R' nK 
(upper panel) and Rx (bottom panel). Symbols as in Fig.@] 

dances derived in very active stars are likely plagued by sys- 
tematic effects (see Morel & Micela 2004 in the case of oxy- 
gen). Until these effects are better understood, great caution 
must therefore be exercised when using these data. 

It is remarkable that all our programme stars (see also Paper 
II) share about the same evolutionary status, i.e. they are at 
the bottom of the red giant branch or are starting to ascend 
it (Fig. |2j. This supports the idea that enhanced chromospheric 
activity is associated with a specific evolutionary stage during 
which the spinning up of the stars, coupled with the formation 
of a deep convective envelope, greatly enhances dynamo action 
(Barrado et al. 1994). We should point out, however, that our 
sample is not free from selection effects: we have preferentially 
selected slow rotators, which are presumably characterized by 
a comparatively low level of coronal activity. 
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Appendix A: EW measurements 

Table |A~TI presents the line list, oscillator strengths and EW 
measurements (only available in electronic form). 

Appendix B: NLTE corrections 

Table IB. II gives the NLTE corrections applied to the abun- 
dances of Mg, Na and Li. For completeness, we also include 
the stars analyzed in Paper II. 
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Table A.l. Line list, calibrated atomic data and EW measurements. The adopted solar abundances are quoted for each element. 
A blank indicates that the EW was not reliably measurable, either because the line was affected by cosmic rays/telluric features 
or significantly van der Walls broadened (see Paper II), or because the Gaussian fit was judged unsatisfactory. 
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33.4 








29.2 


50.7 










Tii; log e e = 


=4.99 




























5766.330 




3.294 


0.370 


37.5 


31.6 


35.0 


35.1 


31.2 


30.3 


61.3 


36.3 




46.7 


44.2 


Cri; log eo- 


=5.67 




























5787.965 




3.323 


-0.138 


81.7 


77.9 


81.2 


81.6 


84.5 


72.7 


108.0 


88.3 




98.9 


97.4 


6882.475 




3.438 


-0.238 


69.3 
















109.9 




: 


6882.996 




3.438 


-0.305 


63.3 
















97.4 






6925.202 




3.450 


-0.227 








78.2 


72.0 




107.3 








94.7 


Fei; log e e = 


=7.67 




























5543.937 




4.218 


-1.155 


88.1 


92.6 


93.9 


92.2 


91.0 


86.7 


110.2 


94.3 


122.6 


100.8 


101.2 


5638.262 




4.221 


-0.882 


106.5 


111.1 


115.5 


113.5 


108.7 


100.7 


134.0 


109.6 


155.4 


119.5 


112.0 


5679.025 




4.652 


-0.863 


79.7 


85.5 


83.2 




78.2 






85.0 




88.9 


87.5 


5732.275 




4.992 


-1.473 


37.1 


43.7 




42.7 




39.5 


52.7 


34.2 


40.3 


38.8 


43.6 


5806.717 




4.608 


-0.984 


83.1 


88.2 


88.0 


91.1 


76.4 


77.3 


105.7 


83.7 




91.1 


90.6 


5848.123 




4.608 


-1.282 




77.7 




77.4 


65.8 


68.3 






103.6 


78.0 


80.3 


5855.091 




4.608 


-1.681 


44.9 


49.3 


48.9 


48.3 




42.2 


65.5 


41.8 


64.0 


47.5 




5905.689 




4.652 


-0.860 


81.2 






92.3 






100.2 






88.8 


88.0 


5909.970 




3.211 


-2.731 








75.9 


82.4 


69.0 




76.4 








5927.786 




4.652 


-1.243 


65.7 


70.5 


69.6 


66.1 


57.7 


63.5 




64.1 


87.7 


68.4 




5929.667 




4.549 


-1.332 


69.5 






69.9 


62.8 


66.3 


86.4 


67.2 




73.2 


82.6 


5930.173 




4.652 


-0.347 


112.7 


116.0 




122.9 


112.3 


106.4 


134.8 


117.0 


160.6 


126.5 


127.1 


5947.503 




4.607 


-2.059 












23.1 


46.1 






28.5 




6078.491 




4.796 


-0.414 


97.5 


104.3 


106.8 


112.3 


93.7 


97.1 


121.9 


104.0 


126.1 


110.8 


110.6 


6078.999 




4.652 


-1.123 


71.8 


73.3 


75.2 


77.1 


66.9 


70.3 




74.0 


95.3 


79.4 




6094.364 




4.652 


-1.749 


43.6 


43.2 


44.7 


50.6 


31.1 


35.7 


61.0 


38.3 




47.6 




6098.280 




4.559 


-1.940 


41.4 




43.2 


43.4 








38.5 




45.6 




6151.617 




2.176 


-3.486 


98.6 


98.3 


99.2 


89.3 


107.1 


85.4 


126.6 


104.0 


149.6 


114.1 


113.7 


6165.361 




4.143 


-1.645 


77.8 


78.1 


78.6 


77.3 


76.8 


71.9 


96.9 


76.2 


119.7 


86.3 


86.8 


6187.987 




3.944 


-1.740 


81.1 


83.8 


89.6 


82.8 


78.5 


76.2 


103.9 


83.3 


119.6 


93.1 


94.3 
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Table A.l. Continued. 



-1(A) A-(eV) log gf EW(mA) 

HD28 HD 1227 HD 4482 HD 17006 HD 19754 HD 154619 HD 156266 HD 181809 HD 182776 HD 202134 HD 204128 HD 20 



6219.279 


2.198 


-2.551 


147.0 


149.9 


156.1 


142.3 


164.9 


132.3 












6252.554 


2.404 


-1.867 












164.0 












6322.690 


2.588 


-2.503 


127.8 


125.1 


129.9 


120.9 


134.3 


112.6 


161.5 


130.1 






146.9 1 


6335.328 


2.198 


-2.432 


159.0 


152.9 


166.6 


151.8 


173.2 


141.7 












6336.823 


3.687 


-0.896 


141.8 


139.4 


150.4 


154.5 


152.6 




169.7 








171.1 


6436.411 


4.187 


-2.538 


36.7 


37.9 




35.2 








27.3 








6593.871 


2.433 


-2.342 










159.6 










173.8 




6699.162 


4.593 


-2.172 


27.6 






31.3 




26.1 


45.5 


22.3 




25.3 




6713.771 


4.796 


-1.606 


46.7 


46.3 


47.4 


52.9 


36.4 




65.5 


41.9 




49.3 




6725.353 


4.104 


-2.370 


43.1 


45.7 


47.9 


44.5 




38.3 


63.4 


39.5 




43.9 




6726.661 


4.607 


-1.200 


70.7 


78.8 


75.5 


76.3 


70.9 


67.5 


91.1 


72.0 




76.1 


80.0 


6733.151 


4.639 


-1.594 


52.6 


58.0 


54.5 


56.2 


47.5 


51.5 


72.8 


48.1 




57.1 


59.3 


6745.090 


4.580 


-2.192 


23.6 




27.6 






21.7 




24.3 








6750.150 


2.424 


-2.727 


124.1 


127.7 


128.0 


120.3 


135.4 


114.6 


160.2 


133.7 




146.5 


151.9 1 


6806.847 


2.728 


-3.265 


84.8 




82.7 


75.4 


85.0 


73.3 




83.5 


129.7 


99.1 




6810.257 


4.607 


-1.129 


73.3 


83.6 


81.3 


85.9 


74.5 


74.8 


98.8 


79.7 


117.8 


84.7 


94.5 1 


6820.369 


4.639 


-1.289 


71.3 


76.1 


76.1 


76.0 


60.4 


68.3 




68.8 




78.7 




6843.648 


4.549 


-0.934 


86.2 


95.3 


94.6 


97.2 


82.8 


86.4 


111.8 


91.4 




94.5 


94.5 1 


6857.243 


4.076 


-2.203 


51.7 


53.2 


53.1 


50.2 




48.7 


73.4 


47.2 






53.8 


6862.492 


4.559 


-1.509 


56.4 


62.8 


61.5 


64.1 


49.2 


62.5 


80.8 






64.7 


64.2 


6988.523 


2.404 


-3.623 












74.0 












7022.953 


4.191 


-1.184 


92.9 




99.0 












133.4 


104.4 


! 


7219.678 


4.076 


-1.715 








78.9 




73.3 


97.5 






84.0 




7306.556 


4.178 


-1.684 








78.9 






100.2 






86.6 


84.8 


7746.587 


5.064 


-1.379 


43.3 


47.5 




50.8 












46.1 




7748.274 


2.949 


-1.748 


163.5 


162.7 


170.2 


161.3 


164.3 


147.1 




168.4 






187.8 1 


7751.137 


4.992 


-0.841 


77.2 


82.7 




82.0 


70.7 


75.7 




71.3 




82.8 


91.9 


7780.552 


4.474 


-0.175 


142.9 


148.4 


152.2 


160.3 


149.1 


140.6 


168.9 


155.8 


193.0 


162.1 


166.4 1 


7802.473 


5.086 


-1.493 


35.9 


37.2 


38.4 


39.4 




33.3 








36.1 




7807.952 


4.992 


-0.602 


88.6 


98.8 


93.2 


95.5 




85.4 


105.5 


86.4 


128.1 


94.8 


103.0 1 


8239.127 


2.424 


-3.433 








91.3 


103.8 


82.5 








117.4 




8592.945 


4.956 


-0.567 








102.9 




93.0 


114.5 










Feii;log<r e =7.67 




























5991.376 


3.153 


-3.702 


46.7 






54.0 




52.7 


58.2 


41.2 


61.2 


46.3 


46.8 


6149.258 


3.889 


-2.858 


38.5 


60.9 


54.7 


50.2 


41.8 


54.0 


53.5 


45.5 




46.2 




6369.462 


2.891 


^.192 


30.6 


41.1 




36.4 




38.3 












6416.919 


3.892 


-2.750 


47.4 


63.8 


59.4 


57.4 


50.2 


58.1 




51.8 


66.1 


54.9 


59.2 


6456.383 


3.904 


-2.209 


64.2 


84.4 


86.7 


78.9 


78.3 


79.4 


76.3 


70.8 


88.3 


75.8 


74.5 


7711.723 


3.904 


-2.625 








66.7 


48.7 


67.4 


62.1 


52.4 




56.5 




Co i; log e =4.92 




























6454.990 


3.632 


-0.233 




46.9 




50.2 


38.3 


37.5 




40.8 




52.7 


56.1 


Nii;log<r e =6-25 




























5593.733 


3.899 


-0.683 


74.3 


71.9 


70.1 


80.5 


59.9 


64.9 


90.7 


65.5 


77.6 


74.4 


81.0 


5805.213 


4.168 


-0.530 


61.9 


68.8 


64.5 


73.3 


55.1 


59.5 


83.9 


61.6 


81.9 


67.8 




6111.066 


4.088 


-0.785 


60.7 


63.7 


61.3 


68.2 


47.4 


53.7 


77.6 


55.6 


80.6 


63.5 


65.0 


6176.807 


4.088 


-0.148 


89.6 


97.7 


94.6 


98.3 


84.6 


85.7 


111.7 


88.1 




98.6 


96.8 1 


6186.709 


4.106 


-0.777 


54.7 


62.9 


63.7 


67.6 




53.2 


76.5 


54.9 




68.5 


62.5 


6204.600 


4.088 


-1.060 


48.9 


52.9 


51.4 


54.9 


36.5 


42.0 


72.9 


39.7 




51.3 


46.0 


6223.981 


4.106 


-0.876 


52.4 


57.2 


58.7 


61.9 




51.8 


78.9 


52.0 




62.9 


59.7 


6772.313 


3.658 


-0.890 


74.4 


81.8 


76.2 


80.8 


75.8 


69.6 


102.7 


77.1 


116.4 


86.8 




7555.598 


3.848 


0.069 


127.1 


132.0 


130.1 


137.9 


120.7 


115.5 


157.8 


126.1 


160.9 


135.7 


140.4 1 


7797.586 


3.899 


-0.144 


106.6 


111.3 


113.9 


112.4 


107.4 


102.3 


125.6 


108.4 




113.6 


118.4 1 


Ban;loge e =2.13 




























5853.668 


0.604 


-0.758 


106.7 


121.2 


126.7 


91.1 


144.6 


116.2 


129.1 


127.4 


162.3 


134.1 


132.2 ; 
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Table B.l. NLTE corrections applied to the abundances of Mg, Na and Li (Gratton et al. 1999; Carlsson et al. 1994). The 
corrections are defined as: Ae = log(e) NLTE - log(e) LTE . 



Star 


Mgi^5711 


Nai^6154 


Li i ,16707 


HD28 


+0.06 


+0.02 


+0.22 


HD 1227 


+0.06 


+0.02 


+0.15 


HD 4482 


+0.06 


+0.00 


+0.17 


HD 10909 


+0.05 


-0.02 


+0.20 


HD 17006 


+0.08 


-0.04 


+0.10 


HD 19754 


+0.03 


-0.09 


+0.24 


HD 72688 


+0.06 


-0.03 


+0.15 


HD 83442 


+0.06 


-0.06 


+0.24 


HD 113816 


+0.06 


-0.03 


+0.24 


HD 118238 




-0.09 


+0.28 


HD 119285 


+0.06 


-0.03 


+0.22 


HD 154619 


+0.10 


-0.09 


+0.12 


HD 156266 


+0.04 


-0.06 


+0.26 


HD 181809 


+0.06 


+0.02 


+0.19 


HD 182776 




-0.09 


+0.25 


HD 202134 


+0.05 


-0.03 


+0.23 


HD 204128 


+0.06 


-0.04 


+0.20 


HD 205249 


+0.06 


-0.06 


+0.13 


HD 211391 


+0.06 


-0.03 


+0.16 


HD 217188 


+0.06 


-0.01 


+0.20 


HD 218527 


+0.06 


+0.03 


+0.15 



